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by 
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SUMMARY 


Aureole  2  rocket  ^197 3— 1 0 7B)  was  launched  on  26  December  1973  into  an 
orbit  of  inclination  74^  and  eccentricity  0. I  and  has  an  estimated  lifetime  of 
30  years.  The  orbit  has  been  determined  from  observations  for  90  epochs  between 
September  1983  and  December  1984,  during  which  time  the  orbit  was  expected  to  be 
influenced  significantly  by  the  effects  of  27:2  resonance  with  the  Earth's  gravi¬ 
tational  fieldj  exact  resonance  occurred  on  28  April  1984.  -The' observations 
numbered  nearly  7400,  of  which  344  were  from  the  Hewitt  cameras  of  the  University 
of  Aston  which  are  sited  at  Herstmonceux  in  England,  and  Siding  Spring  in 
Australia.  The  orbital  inclination  and  eccentricity  of  the  orbits  derived  had 
standard  deviations  corresponding  on  average  to  positional  accuracies  of  1 30  m 
cross-track  and  80  m  in  perigee  distance. 

The  variations  in  inclination  and  eccentricity  have  been  analysed  individu¬ 
ally  to  determine  values  of  two  pairs  of  lumped  harmonics  of  order  27  from  each 
parameter;  when  these  parameters  were  fitted  simultaneously  they  gave  three  pairs 
of  harmonics  with  standard  deviations  corresponding  to  accuracies  of  approximately 
2.5  cm  in  geoid  height. 
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1  INTRODUCTION 

When  the  satellite  Aureole  2  was  launched  on  26  December  1973,  its  rocket, 
designated  1973-107B,  entered  an  orbit  with  an  estimated  lifetime  of  30  years. 

The  rocket  is  cylindrical  in  shape,  7.4  m  long  with  a  diameter  of  2.4  m,  and  has 
a  mass  of  about  2200  kg.  Its  initial  orbital  parameters*  were:  inclination  74.01°, 
perigee  and  apogee  heights  396  and  1965  km  respectively,  and  nodal  period 
109.02  min. 

In  April  1984,  Aureole  2  rocket  passed  through  the  ron-fi  n'np  of  ?7*?  reson¬ 
ance,  ie  the  track  over  the  Earth  repeated  every  2  days  after  27  revolutions.  If 
the  passage  through  resonance  of  an  orbiting  body  is  slow  enough,  the  effects  of 
27th-order  harmonics  in  the  geopotential  can  build  up  and  result  in  an  appreciable 
perturbation  to  some  of  the  orbital  elements.  Thus  measurement  of  these  resultant 
perturbations  provides  a  good  method  for  accurately  determining  the  appropriate 
lumped  geopotential  harmonics.  The  aim  of  this  Report  is  to  compute  accurate 
orbits  from  observations  made  during  the  time  when  the  27:2  resonance  with  the 
geopotential  was  affecting  the  orbit  and  to  evaluate  lumped  geopotential  harmonics 
of  order  27  from  the  changes  they  produced  in  the  orbital  inclination  and  eccen¬ 
tricity;  this  is  the  first  occasion  on  which  27th-order  harmonics  have  been 
obtained  from  resonant  satellite  orbit  analysis.  The  orbit  was  determined 

between  September  1983  and  December  1984  from  radar  and  optical  observations 

2 

using  the  RAE  orbit  refinement  program  PROP,  in  the  PR0P6  version  . 

2  OBSERVATIONS  AND  ORBITS 
2 . 1  Data  sources 

The  orbit  of  1973-107B  has  been  determined  at  90  epochs  between  18  September 
1983  and  22  December  1984  from  7383  observations,  not  including  those  rejected  in 
the  orbit  determinations. 

These  observations  came  from  four  different  sources,  the  most  accurate  being 
those  from  the  University  of  Aston's  Hewitt  cameras  at  the  Royal  Greenwich 
Observatory,  Herstmonceux,  and  at  Siding  Spring  in  Australia;  344  of  these  obser¬ 
vations  were  used  in  27  of  the  90  orbits.  The  second  group  consisted  of  496 
visual  observations  made  by  volunteer  observers  reporting  to  the  Earth  Satellite 
Research  Unit  at  the  University  of  Aston.  The  third  and  largest  group,  of  4200 
radar  observations,  were  made  by  the  US  Navy  Navspasur  system,  kindly  supplied  by 
the  US  Naval  Research  Laboratory  and  the  fourth  group  consisted  of  2343  radar 
observations  from  the  tracking  station  at  RAF  Fylingiiales. 
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2.2  Observational  accuracy 

The  rms  residuals  of  the  observations  have  been  calculated  using  the  RAE 
3  .  . 

computer  program  ORES  ,  and  have  been  distributed  to  the  observers.  Table  1 
gives  the  residuals  for  selected  observing  stations  with  at  least  five  observa¬ 
tions  accepted  in  the  final  orbit  determinations.  The  US  Navy  observations  from 
station  29  are  geocentric,  and  if  they  were  given  in  the  same  form  as  the  topo- 
centric  observations,  their  angular  rms  residuals  would  increase  by  a  factor  of 


Table  1 

Residuals  for  selected  stations 


Station 

Number  of 
accepted 
observations 

Rms  residuals 

Range 

km 

Minutes  of 

arc 

RA 

Dec 

Total 

1  US  Navy 

497 

1  .5 

1  .3 

2.0 

2  US  Navy 

441 

3.3 

3.  1 

4.5 

3  US  Navy 

456 

3.0 

2.3 

3.8 

4  US  Navy 

461 

3.3 

2.4 

4.1 

5  US  Navy 

481 

2.2 

1.9 

2.9 

6  US  Navy 

516 

1.6 

1.5 

2.2 

29  US  Navy 

1348 

0.6 

0.2* 

0.2* 

414  Capetown 

39(37) 

1  .2 

1.4 

1  .8 

21 15  Yateley 

15(13) 

4.3 

3.  1 

5.3 

2122  Malvern  5 

6(6) 

2.4 

2.1 

3.2 

2265  Farnham 

61(55) 

2.5 

2 . 3 

3.4 

2392  Cowbeech 

7(7)  | 

0.6 

1  .  7 

1.8 

2414  Bournemouth 

173(152) 

2.4 

2.6 

3.5 

2418  Sunningdale 

13(10) 

0,9 

2.2 

2.4 

2420  Willowbrae 

82(75)  ! 

2.7 

2.8 

3.9 

2430  Stevenage  4 

13(11)  ; 

0.7 

2.  1 

2 . 2 

2539  Dymchurch 

31(28) 

1  .3 

1  .2 

1  .7 

2657  Bridgwater 

14(13) 

1  .6 

2.0 

2.5 

2659  Herstmonceux  3’^ 

219(203)  | 

0.06 

0.05 

0.08 

4156  Apeldoorn 

5(5) 

1  .7 

2.2 

2.7 

4160  Achel  1 

9(9) 

5.4 

3.7 

6.5 

8517  Sacramento 

15(14) 

2.9 

3.5 

4.6 

9652  Siding  Spring"!" 

... 

125(117) 

0.07 

0.07 

0.  10 

*  Geocentric  t  Hewitt  cameras 


NB  Figures  in  brackets  indicate  the  number  of  observations  used  to  calculate 
the  rms  residuals;  ie  those  observations  with  residuals  less  than  twice  the 
rms  value. 
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between  5  and  10.  In  calculating  the  rms  residuals  for  the  visual  observers, 
observations  with  residuals  greater  than  twice  the  rms  have  been  omitted,  the 
numbers  used  being  shown  in  brackets.  This  gives  a  truer  impression  of  the  nor¬ 
mal  accuracy  of  the  observer,  as  it  eliminates  observations  marred  by  poor  visi¬ 
bility  and  possible  deficiencies  in  orbital  fitting. 

The  rms  residuals  of  the  Hewitt  cameras  are  5  seconds  of  arc  from  203 
observations  by  the  Herstmonceux  camera  and  6  seconds  of  arc  from  I  17  observa¬ 
tions  by  the  Siding  Spring  camera.  Since  the  residuals  combine  the  orbital  and 
observational  errors,  and  the  orbital  model  is  less  accurate  than  ihe  observa¬ 
tions,  the  observational  errors  of  the  Hewitt  cameras  are  likely  to  be  less  th -a 
their  rms  residuals,  and  2  seconds  of  arc  would  be  an  accuracy  consistent  with 
the  results. 

2 . 3  Orbits  and  orbital  accuracy 

Orbits  were  determined  using  RAE’s  orbit  refinement  program  PROP,  in  the 
PR0P6  version,  and  orbital  elements  at  the  90  epochs  together  with  their  stan¬ 
dard  deviations  are  listed  in  Table  2  on  page  16.  The  epoch  for  each  orbit  is 
at  00  hours  on  the  day  indicated,  and  the  PROP  program  fits  the  mean  anomaly  M 
by  a  polynomial  of  the  form 

M  =  Hq  +  Mjt  +  M2t2  +  M3t3  +  M^t4  +  M5t5  ,  (1) 

where  t  is  the  time  measured  from  epoch  and  the  number  of  M-coef ficients  used 
depends  on  the  drag.  For  1973  107B,  with  orbital  eccentricity  approximately  0.1, 
and  perigee  and  apogee  heights  of  about  400  and  1700  km  respectively, 
were  sufficient  for  87  of  the  90  orbits.  The  other  three  orbits  required  the  use 
of  coefficients  . 

The  value  of  e  ,  the  parameter  which  indicates  the  measure  of  fit  of  the 
observations  to  the  orbit,  varied  between  0.33  and  0.83  with  an  average  value  of 
0.56,  showing  that  all  of  the  orbits  were  fitted  satisfactorily.  The  average 
number  of  observations  in  an  orbit  determination  was  82,  spread  over  a  time 
interval  averaging  4.9  days. 

The  average  standard  deviation  in  eccentricity  e  for  the  90  orbits  is 
0.00001 1 ,  equivalent  to  an  error  in  perigee  distance  of  80  m;  the  average  for  the 
27  orbits  containing  Hewitt  camera  observations  was  0.000008.  The  perigee  dis¬ 
tances  Q  =  a(l  -  e)  from  Table  2  are  plotted  in  Fig  1,  and  exhibit  the  usual 
sinusoidal  oscillations  dependent  on  the  argument  of  perigee,  m  ;  also  plotted 

1   
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in  Fig  1,  but  on  a  larger  scale,  are  the  values  of  Q'  ,  the  perigee  distance 
after  removal  of  lunisolar  and  zonal  harmonic  perturbations. 

The  mean  standard  deviation  in  inclination  for  the  90  orbits  is  0.0010°, 
corresponding  to  an  error  of  130  m  in  cross-track  distance;  for  the  27  orbits 
containing  Hewitt  camera  observations  the  accuracy  is  again  better,  the  average 
standard  deviation  being  0.0008°. 


3  THEORY  FOR  THE  RESONANCE  EFFECTS 


The  theory  has  been  g'ven  in  detail  in  Ref  4  and  will  only  be  summarized 
here.  The  longitude-dependent  part  of  the  geopotential  at  an  exterior  point 
(r,0,X)  is  written  as^ 


co  £ 

;  E  m 


(cos 


e) 


C„  cos  mA 
Am 


+  S„  sin  mA>  N. 
Am  I  Jim 


(2) 


A=2  m= 1 


where  r  is  the  distance  from  the  Earth's  centre,  0  is  co-latitude,  A  is 

longitude  (positive  to  the  east),  p  is  the  gravitational  constant  for  the 
3  2 

Earth  (398600  km  /s  ) ,  R  is  the  Earth's  equatorial  radius  (6378.1  km), 

P™  (cos  0)  is  the  associated  Legendre  function  of  order  m  and  degree  A  ,  and 

C.  and  S  are  the  normalized  tesseral  harmonic  coefficients,  of  which  only 
Jim  Jim 

those  of  order  m  =  27  are  relevant  here.  The  normalizing  factor  N^m  is 
given  by 

„2  =  2(2£  +  1 ) ( A  -  a)! 

WAm  (A  +  m) !  *  K 


The  rate  of  change  of  inclination  i  caused  by  a  relevant  pair  of  coef- 

.  .  -  -  4  . 

ficients,  C.  and  S.  ,  near  B:a  resonance  may  be  written  (ignoring  terms 
2  *m  £m 

of  order  e  )  as 


di 

dt 


JC  (£\ 

jin  i  \a  / 


F.  G.  (k  cos  i  -  m)(R  j 
Amp  Apq 


£-m+l (- 


Am 


-  jS  )  exp{j(Y$  -  qw) 


>] 


(4) 


where  F.  is  Allan's  normalized  inclination  function  ,  G„  is  a  function 
Amp  Apq 

of  eccentricity  for  which  explicit  forms  and  a  computer  program  are  given  in 
Ref  4,  denotes  'real  part  of'  and  j  =  /-T  .  The  resonance  angle  4>  is 
defined  by  the  equation 


<t> 


a(ai  +  M)  +  8(ft  -  v) 


(5) 
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where  w  is  the  argument  of  perigee,  M  the  mean  anomaly,  fi  the  right  ascen¬ 
sion  of  the  node  and  v  the  sidereal  angle.  The  indices  y,  q,  k  and  p  in 
equation  (4)  are  integers,  with  y  taking  the  values  1,  2,  3  ...  and  q  the 
values  0,  ±1,  ±2,  ...;  the  equations  linking  £,  m,  k  and  p  are:  m  =  Y  e  ; 
k=ya-q;  2p  =  £  -  k  . 

Here  8  =  27  and  a  =  2  ,  and  as  we  shall  only  consider  the  y  =  1  terms, 
which  are  usually  dominant,  we  have  m  =  27  and  k  =  2  -  q  .  The  values  of  £ 
to  be  taken  must  be  such  that  £  ^  m  and  (£  -  k)  is  even.  The  successive 
coefficients  which  arise  (for  given  y  and  q)  may  be  grouped  into  a  lumped 
harmonic,  written  as 


£  £ 


where  £  increases  in  steps  of  2  from  its  minimum  permissible  value  £^  ,  and 
the  are  constant  coefficients,  with  =  1  .  The  values  of  the  can 

be  obtained  from  equation  (4),  and  R.H.  Gooding  has  written  a  computer  program 
PROF  for  their  evaluation. 

The  rate  of  change  of  eccentricity  produced  by  a  relevant  pair  of 

coefficients  C„  and  S  near  8:a  resonance  may  be  written 
£m  £m  J 


. (7) 

2 

where  terms  of  order  e  have  again  been  ignored. 

As  the  G„  are  of  order  e  ^  ,  it  turns  out  that,  for  orbits  of 
£pq 

eccentricity  less  than  0.1,  the  leading  terms  in  equation  (4)  are  those  with 
q  =  0  and  q  =  ±1  ,  while  the  main  terms  in  (7)  are  those  with  q  =  ±1  .  These 

are  the  only  terms  that  will  be  evaluated  in  the  analysis  of  1973-107B. 

The  explicit  forms  of  equations  (4)  and  (7)  are  given  in  Ref  4  for  the  31:2 

and  29:2  resonances,  but  not  for  the  27:2  resonance.  The  equation  for  di/dt  at 

27:2  resonance,  with  q  =  0  and  q  =  ±1  ,  is 


8 


i  m27 

/  R  \  —  ,  1 

[_0, 2 

Csl 

O 

1 

U) 

U  /  F28,27, 13c28, 13,0(27  cosec  1 ~  2  cot 

l'27 

sin  $  +  cos 

+  F27  27  13G27  13  1  ^27  cosec  i_  coti)<C2?  sin(<J>  -  tu)  -  cos(*  -  w) 


+  F27, 27, 12G27, 12,-1 <'27  cosec  1  3  cot 


2,-J 


1,1 

7 

_-'.3 


.U 
11 

_-l,3 


C27  sin  ($>  4  uj)  ~S27  cos(<J>  4  w)> 


(8) 


A  factor  (1  -  e  )  should  be  introduced  on  the  right-hand  side  if  terms  of  order 
2 

e  are  required.  The  equation  for  de/dt  ,  with  q  =  ±1  terms  only,  is 


de 

dt 


f  (lf[- 


.1  ’  1 

'27 

-1,3 


F27 ,27,1 3G27 , 13,1  |C27  Sin($  ~  S27  C°s($  "  w)j 

-1,3 


+  F27,27, 12G27, 12,-1  C27  sin(<1>  +  u)  "  S27  COs(<J>  +  u) 


.  (9) 


Here,  with  ya  =  2  and  q  ^  0  ,  the  factors  that  need  to  be  introduced  on  the 

2 

right-hand  side  of  (9)  to  take  account  of  terms  of  order  e  are  given  by 
2  A 

{1  -  e  (1  4  1/q)  4  0(e  )}  ,  from  Ref  4.  Thus  for  the  first  term  in  curly 

2  ^ 

brackets  in  (9)  (q  =  1),  the  factor  is  { 1  -  2e  4  0(e  )}  ,  and  for  the  second 
term  (q  =  -1)  ,  the  factor  is  { 1  +  OCe^)}  . 


4  ANALYSIS  OF  THE  VARIATIONS  IN  INCLINATION  AND  ECCENTRICITY 

4 . 1  Progress  through  resonance 

Variations  in  $  during  the  period  covered  by  this  study  are  shown  in 
Fig  2.  The  increase  in  <t>  ,  from  about  -10  deg/day  through  to  +7  deg/day,  pro¬ 
ceeds  quite  steadily,  but  is  rather  slower  after  the  exact  resonance,  which 
occurs  on  28  April  1984.  Progress  through  resonance  is  not  as  slow  as  would 
ideally  be  hoped  for:  it  is  8  times  faster  than  for  the  recent  analysis7  of 
1968-40B  at  29:2  resonance,  so  the  resulting  perturbations  and  derived  coef¬ 
ficients  cannot  be  expected  to  be  so  accurate  as  those  obtained  from  1968-40B. 

4.2  Analysis  of  inclination 

Jhe  raw  values  of  inclination  given  in  Table  2  need  to  be  cleared  of  per¬ 
turbations  due  to  zonal  harmonics  and  lunisolar  effects:  this  has  been  done  by 

g 

use  of  the  PROD  computer  program  with  integration  at  1-day  ^“'rvals. 

Perturbations  due  to  the  J  „  harmonic  are  recorded  within  each  PROP  run 
and  have  also  been  removed.  Fig  3  shows  the  resulting  values  of  i  ,  with  sd. 
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These  values  of  inclination  were  then  fitted  with  the  computer  program 
9 

THROE  ,  within  which  the  effects  of  atmospheric  rotation  are  also  removed  (the 
value  used  for  the  atmospheric  rotation  rate  A  was  1.00,  in  conformity  with 
Ref  10).  At  first  the  fittings  were  made  with  all  three  pairs  of  values  of  (y,q) 
in  equation  (4),  that  is,  (y,q)  =  (1,0),  (1,1)  and  (1,-1),  as  in  equation  (8). 
However,  it  was  found  that  the  values  for  (y,q)  =  (1,-1)  were  small  and  indeter¬ 
minate,  so  these  terms  were  dropped.  It  also  became  apparent  from  the  THROE  runs 
that  a  number  of  the  orbits  fitted  badly,  as  indeed  is  obvious  from  Fig  3.  Values 
of  inclination  which  had  residuals  greater  than  2c  ,  where  e  is  the  overall 
measure  of  fit,  were  relaxed  by  doubling  the  standard  deviation  and,  if  necessary, 
quadrupling  it.  Also  the  first  four  orbits  were  omitted.  As  a  result  of  these 
procedures  12  values  had  their  standard  deviations  doubled  and  two  values  had 
quadrupled  standard  deviations.  The  fitting  began  at  MJD  45628,  covered  85  orbits 
and  ended  at  MJD  46051,  the  last  orbit  also  being  omitted,  as  it  was  ill-fitting. 
The  overall  measure  of  fit,  e  ,  was  1.17,  which  is  quite  satisfactory. 

Although  the  simultaneous  fitting  of  i  and  e  will  be  preferred  (see 
section  4.4),  the  values  of  the  (y,q)  =  (1,0)  terms  from  the  fitting  of  i  alone 
should  be  fairly  satisfactory.  They  are: 

0,2  0,2 

10  C  =  22.4  ±  12.4  ,  10  S  =  30’8  1  l0*2  •  (10) 

The  'hats'  over  the  lumped  harmonics  indicate  that  these  are  the  values  emerging 
from  THROE,  within  which  is  replaced  by  an  approximation  *  See 

Ref  4  and  section  4.4.  (The  values  of  (C.S)^'  are  not  giver,  re  they  are  sub¬ 
sidiary  terms  and  are  normally  better  determined  from  e  .) 

4 . 3  Analysis  of  eccentricity 

As  the  eccentricity  of  1973—10 7B  is  quite  large,  and  decreases  appreciably 
due  to  drag  during  the  resonance,  it  is  better  to  work  with  the  perigee  distance, 
which  is  much  less  affected  by  drag.  The  lower  graph  in  Fig  1  gives  the  raw 
values  of  a(l  -  e)  ,  showing  the  characteristic  oscillation  due  to  odd  zonal 
harmonics,  which  has  an  amplitude  near  7  km.  The  upper  graph  in  Fig  1  shows  the 
values  after  removal  of  zonal  harmonic  and  lunisolar  perturbations.  These  values, 
denoted  by  Q'  ,  should  show  the  effects  of  air  drag  and  of  resonance  only.  For 
the  fitting  by  THROE,  it  is  convenient  to  define  a  revised  value  of  e  ,  e 
say,  based  on  Q'  .  Thus 


e 

rev 


-  Q'/a 


(ID 
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where  a  is  the  mean  value  of  a  during  the  orbit  determinations. 

To  a  first  approximation  the  decrease  in  Q'  due  to  drag  at  any  stage  in  the 
orbit  determinations  is  given  by  sHinCe^/e)  ,  where  H  is  the  density  scale 
height,  e^  the  initial  eccentricity  and  e  its  current  value  (corrected  to 
to  =  Wq)  .  On  comparing  values  of  e  at  similar  values  of  a1  in  Table  1,  the 
decrease  in  e  between  MJD  45595  and  45983  is  found  to  be  by  a  factor  1.0112, 
and,  on  taking  H  =  50  km,  this  gives  a  decrease  of  0.28  km  in  Q'  ,  which  trans¬ 
lates  into  an  average  increase  of  almost  exactly  10  7  per  day  in  e  .  Thus  a  new 
value  of  e  ,  e*  say,  cleared  of  air  drag,  may  be  defined  as 

e*  =  e  -  10~7(t  -  t.)  .  (12) 

rev  0 

(In  practice  the  correction  was  taken  as  -|N  *  10  ^  ,  where  N  is  the  orbit 
number . ) 

The  resulting  values  of  e*  were  fitted  by  THROE  with  (v,q)  =  (1,1)  and 
(1,-1).  As  with  i  ,  the  first  four  values  were  omitted,  and  a  number  of  the 
standard  deviations  were  relaxed.  The  two  anomalous  values  at  MJD  46033  and 
46037  (see  Fig  1)  were  allocated  standard  deviations  of  0.0001,  and  eight  other 
standard  deviations  which  exceeded  2e  were  relaxed  by  a  factor  of  2.  The  over¬ 
all  measure  of  fit,  e  ,  then  had  the  value  2.66:  this  is  rather  high,  but  not 
unusually  so,  because  values  of  e  near  2  often  arise  when  fitting  eccentricity7. 

4 . 4  Fitting  of  i  and  e  simultaneously 

As  both  inclination  and  eccentricity  required  the  ,(v,q)  =  (1,1)  term,  it 
seemed  best  to  fit  them  simultaneously  with  R.H.  Gooding's  SIMRES  program.  The 
SIMRES  fitting  was  made  with  (y,q)  =  (1,0), (1,1)  and  (1,-1)  and  the  eccentricity 
fitting  was  given  a  lower  weighting,  in  accordance  with  the  ratio  of  the  values 
of  e  in  the  THROE  fittings  for  i  and  e  :  the  weighting  factor  was  1.848. 

The  overall  measure  of  fit  for  the  SIMRES  fitting  was  1.362,  and  the  individual 
values  of  e  for  i  and  e  were  1.32  and  2.96,  as  compared  with  1.17  and  2.66 
for  the  individual  fittings.  Thus  the  fittings  are  not  much  worse,  and  the  com¬ 
bined  fitting  is  to  be  recommended  because  the  (y,q)  =  (1,1)  term  is  significant 
for  both  i  and  e  ,  and  SIMRES  provides  an  'average'.  The  values  of  the  har¬ 
monics  that  emerge  are: 
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9- 
10  C 


0,2 

27 


9 

10C 


,1 

27 


9*-l,3 
10  C27  = 


g.0,2 

18.4 

+ 

Ui 

00 

10  S27 

q.1,1 

14.4 

± 

15.1  , 

10  S„_, 

27 

q.-l, 3 

14.5 

+ 

8.7  , 

10  S27 

26.3  '  12.3 

-9.1  .+  13.3 

2.4  ±  11.2 


(13) 


-  -  o  2 

It  may  be  noted  tl.at  the  values  of  (C,S)27  in  equations  (10)  are  consistent 
wi th  those  in  (13). 

The  fittings  of  inclination  and  eccentricity  are  shown  in  Figs  4  and  5, 
where  the  standard  deviations  indicated  are  those  after  relaxation.  It  will  be 
seen  that  in  the  combined  fitting  some  of  the  residuals  exceed  2e  :  for  example 
the  residual  for  i  at  MJD  45734  is  3.1(=  2.3e)  .  Further  readjustment  of  the 
relaxations  was  not  attempted. 

The  values  of  C  and  S  above  have  been  given  'hats'  (")  to  indicate  that 

they  are  the  values  emerging  from  SIMRES,  in  which  the  values  of  the  C  functions 

in  equations  (8)  and  (9)  are  replaced  by  an  approximation  G  .  Thus  if  C  is 

m 

the  correct  value  we  have 


.q,k 


G„  C 

Vo’  m 


q,k 


CVo‘>C" 


(14) 


Thus  the  values  of  C  and  S  in  equations  (13)  have  to  be  divided  by 

C-  /G,  ,  values  of  which  (always  >1)  have  been  obtained  from  the  computer 

Voq  £opoq 

program  GQUAD^ .  For  an  orbit  of  such  high  eccentricity  as  1973-107B  and  such 
high  degree,  the  corrections  are  large:  the  three  coefficients  arising  need  to  be 
divided  by  3.323,  2.014  and  2.009  respectively. 


Thus  the  values  given  in  equations  (13)  may  be  rewritten  as 


q  °»2 

10  C  =  5.5  ±  4.8  , 


q_°>2 

10  S2y  =  7.9  ±  3.7 


,9- 


1,1 


I0'C27  =  7.1  ±  7.5  , 

o_-l>3 

10  C2?  =  7.2  ±  4.3  , 


9- 

10  S 


9- 

10  S 


,1 

27 

-1,3 

27 


=  -4.5  ±  6.6 


1.2  ±  5.6 


.  (15) 


These  are  the  final  values  of  the  lumped  harmonics  derived  from  analysis  of 

1973-107B. 
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For  eccentricity,  further  THROE  runs  were  made  in  which  the  value  of  the 
third  zonal  harmonic  was  adjusted  to  minimize  the  value  of  e  ,  and  e  was 

substantially  reduced,  from  2.66  to  1.85.  However,  this  ’optimum-J^’  run  for  e 
and  the  THROE  run  for  i  ,  when  combined  in  SIMRES,  gave  larger  standard  devia¬ 
tions  for  the  lumped  harmonics  than  those  in  (15):  so  the  values  (15)  are 
preferred. 


LUMPED  HARMONICS  IN  TERMS  OF  INDIVIDUAL  COEFFICIENTS 


^Hm’ ^£m 


The  lumped  harmonics  (C,S) 


r 


27 


are  expressible  in  terms  of  the  individual 


coefficients  (C.S)^  by  equations  (6),  and  the  computer  program  PROF  evaluates 
the  Q  functions  with  adequate  accuracy  if  e  is  very  small.  For  1 973-1 07B, 

w  •  •  ok 

however,  e  —  0.087  and  it  is  necessary  to  multiply  each  Qj”  by  a  correction 


factor  £  »  say,  where 


1 1 


£  = 


J£pq  V0q 


Lpq 


=  f 


Voq 


J£pq 

’*pq 


(16) 


and  f  is  the  factor  by  which  the  lumped  harmonics  had  to  be  multiplied,  namely 
0.3009,  0.4965  and  0.4978  respectively,  for  the  three  pairs  of  harmonics  in 
equations  (15).  Values  of  ^£pq^£pq  have  been  obtained  from  the  computer  pro¬ 
gram  GQUAD  for  values  of  £  up  to  48  with  e  =  0.087  and  the  resulting  values 
of  £  are  given  in  Table  3.  It  is  apparent  that  £  departs  greatly  from  I  for 


Table  3 

Values  of  £  for  1973-1 07B  with  e  =  0.087 


q  =  0 

q  =  ±1 

£ 

K 

£ 

£ 

28 

1  .000 

27 

1  .000 

30 

1  .152 

29 

1  .098 

32 

1.331 

31 

1.212 

34 

1.543 

33 

1.341 

36 

1.792 

35 

1.490 

38 

2.085 

37 

1  .662 

40 

2.432 

39 

1  .858 

42 

2.840 

41 

2.084 

44 

3.321 

43 

2.344 

46 

3.889 

45 

2.642 

48 

4.560 

47 

2.986 
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high  values  of  £  .  The  resulting  expressions  for  the  lumped  harmonics  in  terms 


of  the 

individual 

coefficients , 

after  the  values 

of  Q  from 

PROF  have  been 

multiplied  by  £  , 

,  are  as  follows: 

0,2 

C27 

C28,27 

+  0-,39S),27 

-  °'3I9'C32,27  - 

’  °-395£34,27 

-  °-256S6,37 

-  °-06'S8,27  * 

°-089540,27 

-  0.156C42j27 

+  °-U4544,27 

9 

(17) 

M 

C27 

C27,27 

-  I.411C29f27 

-  °-785Sl,27  * 

°-094S3,27 

+  °-541S5,27 

*  0'535E37,27  * 

0.278^39,27 

-  0.013B4]>27 

-  0.204E43i27  - 

0 . 255645  ^  27 

> 

(18) 

_-l, 

C27 

3 

^27,27 

~  0-442^29,27 

-  °-695S,,27  - 

0.508033,27 

-  0.188C35i27 

+  °-09°S7.27  * 

°-244S9,27 

+  0.268C. ,  97 

4 1  ,  z7 

+  °-l9^43,27 

-  O.OS3C45j27 

• 

(19) 

Similar  equations  apply  for  the  S  coefficients.  Equation  (17)  has  been  termi¬ 
nated  after  9  terms  at  £  =  44  ,  after  which  no  numerical  coefficient  exceeds  0.1 
Ten  terms  have  been  included  for  equations  (18)  and  (19):  in  the  neglected  terms 
(£  >45),  no  numerical  coefficient  exceeds  0.2. 

6  DISCUSSION 

6 . 1  Comparison  with  comprehensive  gravity-field  models 

A  number  of  recent  comprehensive  models  of  the  gravity  field  give  values  of 

-■  00  C„_  „,... C„,  which  can  be  substituted  into  equation  (17)  to  evaluate 

)  Z  /  3(J )  /  /  Jo  y  Z  / 

the  lumped  harmonic  00,2  ,  and  similarly  for  S  .  The  models  chosen  for  compari 

27  .711 

son  are,  as  in  previous  such  comparisons  ’  ,  the  Goddard  Earth  Model  10B 

(GEM  10B,  Ref  12),  the  1981  model  of  Rapp'3,  and  the  European  GRIM3-L1  (Ref  14). 
The  newer  models  GEM  T1  and  GEM  T2  are  not  included,  because  the  values  for 
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i  =  28  to  36  are  believed  to  be  less  reliable  than  for  the  other  three  models, 
being  based  on  satellite  data  only.  The  values  of  the  lumped  harmonics  obtained 
are  given  in  Table  4. 

Table  4 

Values  of  lumped  harmonics  from  1973-1Q7B  and  comprehensive  models 


q_0,2 

,0C27 

o_0,2 

10  S27 

Q—  1,1 

10  C27 

o-'.l 

10  S27 

Q  1  ,3 

1°  c27 

G--1, 3 

10  S27 

1973-107B 

5.5  ±  4.8 

7.9  ±  3.7 

7.1  ±  7.5 

-4.5  ±  6.6 

7.2  ±  4.3 

1.2  ±  5.6 

GEM  10B 

-13.4 

3.1 

9.3 

-1.0 

6.5 

0.6 

Rapp  1981 

-5.6 

6.2 

17.6 

-19.5 

10.4 

-8.7 

GRIM  3-L1 

-8.6 

4.4 

28.7 

1 

00 

15.6 

-3.6 

The  estimated  standard  deviations  of  the  values  from  the  models  in  Table  4  are 

mostly  near  ±4  ,  slightly  lower  than  those  from  the  resonance.  The  values  in 

Table  4  vary  rather  widely,  but  over  50%  differ  by  less  than  the  sum  of  the 

standard  deviations,  and  the  resonance  values  and  GEM 1  OB  agree  particularly  well, 

-0  2 

apart  from  the  first  coefficient  (C 2’  ).  The  values  from  Rapp  1981  and 
GRIM3-L1  are  very  similar  to  each  other  (probably  due  to  using  similar  terres¬ 
trial  gravity  data  in  the  solution),  and  in  general  do  not  agree  well  with  the 

-0,2  .  . 

resonance  values,  except  for  S^-j  ,  though  all  agree  to  within  about  twice  the 
sum  of  the  standard  deviations.  Thus  it  can  be  said  that  the  values  from 
1973-107B  are  broadly  consistent  with  the  models,  but  are  not  accurate  enough  to 
provide  any  significant  improvement  on  the  models  as  a  whole. 


6 . 2  Geoid  height  accuracy 

The  error  in  geoid  height  implied  by  the  standard  deviations  a  of  the 

lumped  harmonics  may  be  roughly  estimated  as  Ro/Q  ,  where  R  is  the  Earth's 
^  k  2  2  2  j 

radius  and  Q  =  £,,/£  )  }  ,  the  summation  running  from  up  to  the 

u  _0,2  _0,2  u 

maximum  i  considered  (44  or  45).  For  C„  and  S  the  value  of  Q  is 

—9 

1.09  and,  with  a  =  4.3  *  10  as  the  average,  the  error  in  geoid  height  is 

-1,1  -1,1 

about  2.5  cm.  For  and  the  value  of  Q  is  1.75  and,  with 

o  =  7.1  x  10  9  as  the  average,  the  error  in  geoid  height  is  about  2.6  cm.  For 

— —  1 » 3 


'‘27 


and  S 


-1,3 

27 


-  .  -9 

the  value  of  Q  is  1.26  and,  with  a  =  5.0  x  10  as  the 


average,  the  error  in  geoid  height  is  about  2.5  cm. 
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6 . 3  General  discussion 

This  analysis  of  1973-107B  is  the  first  known  attempt  at  determining  lumped 
harmonics  of  order  27  from  an  orbit  passing  through  27:2  resonance.  However  the 
satellite  was  not  ideal  for  the  purpose,  because  perigee  height  was  below  400  km 
and  hence  drag  effects  were  strong  enough  to  carry  the  orbit  through  resonance 
rather  rapidly  -  about  8  times  faster  than  in  Walker's  recent  analysis  of 
1 968-40B  at  29:2  resonance^.  The  accuracy  achieved  here  was  therefore  expected 
to  be  considerably  poorer  than  in  Walker's  analysis,  especially  as  the  effects  of 
the  resonance  on  an  orbit  are  smaller  when  the  satellite's  altitude  is  greater. 
This  expectation  is  confirmed  by  the  results:  in  Walker's  analysis  the  best  geoid 
height  accuracy  was  0.5  cm  (though  the  errors  in  some  harmonics  were  much  larger) 
here  the  values  are  all  near  2.5  cm. 

The  high  eccentricity  of  the  orbit  also  causes  problems,  and  it  is  possible 
that  the  q  =  ±2  terms  in  equations  (4)  and  (7)  may  be  significant. 

The  results  therefore  show  that  analysis  of  27:2  resonance  is  feasible,  but 
that  a  less  eccentric  orbit  of  considerably  lower  drag  is  needed  to  obtain  values 
of  the  lumped  harmonics  tnat  are  much  better  than  those  available  from  the  com¬ 
prehensive  models  of  the  gravity  field. 

7  CONCLUSIONS 

The  orbit  of  1973-I07B  has  been  analysed  at  90  epochs  from  nearly  7400 
observations,  as  it  passed  through  27:2  resonance  between  September  1983  and 
December  1984.  The  orbital  accuracy  was  good,  corresponding  to  positional 
accuracies  of  1 30  m  cross-track  and  80  m  in  perigee  distance. 

The  variations  of  inclination  and  eccentricity  have  been  analysed  to  deter¬ 
mine  three  pairs  of  lumped  harmonics  of  order  27,  with  accuracies  equivalent  to 
approximately  2.5  cm  in  geoid  height.  These  accuracies  are  of  the  same  order  as 
those  of  values  obtainable  from  comprehensive  gravity  field  models:  to  improve  on 
the  latter  values  calls  for  an  orbit  of  considerably  lower  drag  than  1973-107B, 
i-e  having  a  perigee  height  well  above  400  km. 
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